A post-Keplerian parameter to test gravito-magnetic effects in binary pulsar systems 
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We study the pulsar timing, focusing on the time delay induced by the gravitational field of 
the binary systems. In particular, we study the gravito-magnetic correction to the Shapiro time 
delay in terms of Keplerian and post-Keplerian parameters, and we introduce a new post-Keplerian 
parameter which is related to the intrinsic angular momentum of the stars. Furthermore, we evaluate 
the magnitude of these effects for the binary pulsar systems known so far. The expected magnitude 
is indeed small, but the effect is important per se. 
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PACS numbers: 4.20.-q 95.30.Sf 



I. INTRODUCTION 

The first binary radio pulsar (PSR B 1913+16]was dis- 
covered some 30 years ago by Taylor and Hulse[l(. Since 
those years, a great amount of work has been done, and 
pulsars in binary systems have proved to be celestial labo- 
ratories for testing the relativistic theories of gravity (see 
and references therein) . Indeed, up to the present day, 
Einstein's theory of gravity, General Relativity (GR) , has 
passed all observational tests with excellent results. How- 
ever, even if the aim of the experimental relativists is al- 
ways to achieve a greater precision, we must not forget 
that most of the tests of GR come from solar system ex- 
periments, where the gravitational field is in the "weak" 
regime. On the other hand, it is expected that deviations 
from GR can occur for the first time in the " strong" field 
regime: hence, the solar system experiments are inade- 
quate to this end. On the contrary, the strong gravita- 
tional field is best tested by means of pulsars. Pulsars, 
which arc highly magnetized rotating neutron stars, are 
important both for testing relativistic theories of gravity 
and for studying the interstellar medium, stars, binary 
systems and their evolution, plasma physics in extreme 
conditions. As for the tests of gravity, the recent dis- 
covery of the first double pulsar PSR J0737-3039 0,0 
provided an astonishing quantity of data, which make 
this system a rare relativistic laboratory |5| ■ 

As we pointed out elsewhere , this system in particu- 
lar, and binary pulsars systems in general, could be useful 
for testing the so called gravito-magnetic effects. These 
effects are originated by the rotation of the sources of 
the gravitational field, which gives rise to the presence of 
off-diagonal goi terms in the metric tensor. The gravita- 
tional coupling with the angular momentum of the source 
is indeed much weaker than the coupling with mass alone, 
the so called gravito-electric interaction. In fact, the ra- 
tio between the former and the latter can be estimated 
to be of the order of 
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where Rs — 2GM/c 2 is the Schwarzschild radius of the 
source, M being its mass, and J (the absolute value 



of) its angular momentum. At the surface of the Sun, 
which is the most favorable place in the solar system, 
evaluation of eq. (JTJ gives e ~ 10~ 12 , thus evidencing 
the weakness of the gravito-magnetic versus the gravito- 
electric interaction. The smallness of e is the reason 
why, though having been suggested from the very begin- 
ning of the relativistic age [7j , the experimental verifica- 
tion of the existence of gravito-magnetic effects has been 
very difficult until today (see and references therein). 
The relevance of pulsar systems for the detection of the 
gravito-magnetic effects lays in the fact that the ratio (JTJ 
can be less unfavorable whenever r is approaching the 
Schwarzschild radius of the source: this can be the case 
of a source of electromagnetic (e.m.) signals orbiting 
around a compact, collapsed object. 



As for the rotation effects in pulsars binary systems, 
in general, the coupling of the intrinsic angular momen- 
tum of the stars with the orbital angular momentum, 
and the coupling of the intrinsic angular momenta of the 
two stars themselves were studied, together with the cor- 
responding precession effects 0];E3'E3iE3iE3- The 
gravitomagnetic effects on the propagation of light in a 
binary system were studied by Kopeikin and Mashhoon 
(see [l4| and references therein). 



In a previous paper |fj] we studied the effects of the 
gravitational field on e.m. pulses propagating in a bi- 
nary system, and we emphasized the gravito-magnetic 
contribution. In doing that, we used a simplified model 
which considered circular orbits only. Here wc generalize 
those results, by taking into account orbits with arbi- 
trary eccentricity, which is a more realistic approach on 
the basis of the knowledge of binary systems discovered 
until today. In particular, we study the gravito-magnetic 
correction to the Shapiro contribution to the time of flight 
of the signals, focusing on its effect on the arrival times 
perceived by the experimenter on the Earth, and we in- 
troduce a new Post-Keplerian parameter which is related 
to the intrinsic angular momentum of the stars. Finally, 
we evaluate the magnitude of these effects for the binary 
pulsar systems known so far. 
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II. THE CONTEXT: PULSARS TIMING 

When studying pulsars, what is measured are the pulse 
arrival times at the (radio) telescope over a suitably long 
period of time. In fact, even though individual pulses 
are generally weak and have an irregular profile, a reg- 
ular mean profile is obtained by averaging the received 
pulses over a long time. Pulses traveling to the Earth 
are delayed because of the dispersion caused by the in- 
terstellar medium. Besides this delay, other gravitational 
factors influence the arrival times of the signals emitted 
by a pulsar in a binary system: the strong field in the 
vicinity of the pulsar, the relatively weak field between 
the two compact objects forming the binary system and, 
finally, the weak field of the Solar System. Consequently, 
the arrival time Tjy of the iVth pulse, as measured on 
the Earth, depends on a set of parameters ol\,o.2-;Olk, 
which include a description of the orbit of the binary 
system: 

T N = F(N,a u a 2 ,...,a K ) (2) 

In particular, the set a.\,a.2---, <%k includes the Keplerian 
parameters, together with the so called "post-Keplerian" 
(PK) parameters which describe the relativistic correc- 
tions to the Keplerian orbit of the system. Since seven 
parameters are needed to completely describe the dynam- 
ics of the binary system (see [l5j and references therein) , 
the measurement of any two PK parameters, besides 
the five Keplerian ones, allows to predict the remain- 
ing PK parameters. For instance, if the two masses are 
the only free parameters, the measurement of three or 
more PK parameters over-constrains the system and in- 
troduces theory-dependent lines in a mass-mass diagram 
that should intersect, in principle, in a single point |16| . 
This is of course true as far as the intrinsic angular mo- 
menta are not taken into account. 

It is possible to obtain a relation which links the time 
of arrival of a pulse on the Earth to its time of emission. 
More in detail, the following timing formula holds, which 
relates the reception (topocentric) time TEarth on the 
Earth with the emission time T pu i sar in the comoving 
pulsar frame |l7|: 

Tpuisar —TEarth — h — -j^ + Afl + A# Q — A_s Q 

- (Ajj + A E + A s + A a) (3) 

where to is a reference epoch, D/f 2 is the dispersive 
delay (as a function of the frequency of the pulses, /), 
Ar, Ae, As are, respectively the Roemer delay, the Ein- 
stein delay and the Shapiro delay due to the gravitational 
field of the binary system (while A^ ,Ae q , Ag are the 
corresponding terms due to the Solar System field) and 
Aa is the delay due to aberration. 

Here we are concerned with the gravito-magnetic cor- 
rections due to the intrinsic angular momentum of the 
stars to the Shapiro delay Ag, which we analyze in the 
following section. 
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FIG. 1: Notation used for describing the pulsar orbit. 

III. THE SHAPIRO TIME DELAY 

We want to calculate the relation between the coordi- 
nate emission time t e and the coordinate arrival time t a 
(as measured at the Solar System center of mass). 

To this end, we consider a reference frame at rest in 
the center of mass of the binary system. In this reference 
frame, the vector pointing to the pulsar emitting e.m. 
signals is a?i, while the one pointing to its companion star 
is X2', in the following, the suffix "1" will always refer to 
the pulsar, and "2" to its companion. Furthermore, Xf,, 
is the position of the center of mass of the solar system. 

We use the notation of Figure for the description 
of the pulsar orbit. We choose a first set of Cartesian 
coordinates {x,y,z}, with origin in the center of mass 
of the binary system, and such that the line of sight is 
parallel to the z axis. Then, we introduce another set of 
Cartesian coordinates {X,Y, Z}, with the same origin: 
the X axis is directed along the ascending node, the Z 
axis is perpendicular to the orbital plane. The angle 
between the x and X axes is ft, the longitude of the 
ascending node, while the angle between the z and Z axes 
is i, the inclination of the orbital plane. Let x\ = r\X\ 
be the orbit of the pulsar: it is described by 

X\ = cos(w + (p)X + s'm(uj + ip)Y, (4) 

in terms of the argument of the periastron, to, and the 
true anomaly, ip. Let us pose 

9 = to + tp, (5) 
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for the sake of simplicity. Then, we use the notation 

r = x 1 -x 2 (6) 

to describe the position of the pulsar with respect to its 
companion, and we remember that we have, for the Ke- 
plerian problem 



a(l-e 2 ) 



(7) 



1 + e cos ip ' 

where a is the semi-major axis of the relative motion 
and e is the eccentricity. The astronomical elements f2, 
i, uj, a, e represent the Keplerian parameters. In what 
follows, we will use also the definitions r& = |a;&|, r = |r| 
and n = 

rp 

That being said, let us focus on the physical situation 
we are dealing with. If the gravito-magnetic effects are 
neglected, the metric describing the gravitational field of 
the binary system is given by (see, for instance, [Isj ) 

ds 2 = (1 + 2cj))dt 2 - (1 -20) \dx\ 2 , (8) 

where 

Mi M 2 



(x,y,z) = —. 



\x — Xi\ 



\X — X2\ 



(9) 



We see that the total gravitational potential is the sum 
of the two contributions, due to the pulsar, whose mass is 
Mi, and to its companion, whose mass is Afo|2g|. Start- 
ing from © , the mass contribution to the time delay can 
be evaluated following the standard approach, described, 
for instance, in |17| . However, we want to generalize this 
approach in order to take into account the effects of the 
rotation of the sources of the gravitational field, i.e. the 
gravito-magnetic effects. 

To this end, we may guess that the metric, in the co- 
ordinates {X, Y", Z}, assumes the form 

ds 2 = (1 + 2(j>) dt 2 - (1 - 2(f)) (dX 2 + dY 2 + dZ 2 ) 

+ AA ■ dXdt, (10) 

where A is the (total) gravito-magnetic vector potential 
of the system. Eq. (|1(J|I generalizes the weak field metric 
around a rotating source (see ^3j|13)- We suppose that 
the dominant contribution to the total gravito-magnetic 
potential is due to the intrinsic angular momenta Ji , J 2 
of the stars; we suppose also that Ji 7 J 2 are aligned with 
the total orbital angular momentum L , i.e. perpendic- 
ular to the orbital plane. Furthermore, since we assume 
that the signals emitted by pulsar 1 propagate along a 
straight line, and because the gravito-magnetic coupling 
depends on the impact parameter, we may conclude that 
the only relevant gravito-magnetic contribution comes 
from the companion star 2. Consequently, the metric 
(I10II becomes 

ds 2 = (1 + 2<j)) dt 2 - (1 - 20) (dX 2 + dY 2 + dZ 2 ) 

+ ij 2 ( £ - s *}\*dYdt u 2 ( 4^MSdxdt. 



x 2 \ 



x 2 \ 



(11) 



In (II 1|) | a: | varies along the straight line defined by 
x(t) = xi(t e ) + — ^ y- (x b (t a ) - xi(t e )) . (12) 



If we set 



t-t e 
t -t ' 



(13) 



then a — when t = t e , and a = 1 when t = t a . Conse- 
quently, we may write 

x(t) = xi(t e ) + a (x b {t a ) - xi{t e )) . (14) 

After some straightforward manipulations, the line ele- 
ment (|ll|l can be written in the form 



ds 2 = g tt dt 2 + g aa da 2 + 2g ta dtda, (15) 



where 



9tt =1 



2Mo 



Jaa — I 1 



9ta — — 2J 2 



2Mi 

X - Xi\ \x - X 2 \' 

2Mi 2M 2 



\x — Xi\ \x — x 2 \ 



rrb sin i cos ( 



( r b + r i + 2ri rb sin i sin 9) , 
(16) 



\x - x 2 \ 



By setting ds 2 — 0, we easily see that the propagation 
time is made of three contributions, up to first order in 
the masses of the stars, and in the spin angular momen- 
tum: 



At = A< + Ai M + At j, 



(17) 



where 



At Q = [ \Jrl + r\+ 2r b ri sin i sin 6da, (18) 



At 



2M 2 



\x - x 2 \ 



*6 + r i + ^ r i r b sinisin^da, 
(19) 



At; 



rrb sin I cos a 
2J 2 — j-5 —^—da. 

\x - x 2 r 



(20) 



Let us comment on l|18ll - H20f) : the first contribution is a 
purely geometric one, and is due to the propagation, in 
flat space, of e.m. signals from the pulsar toward the So- 
lar System. The second contribution represents the time 
depending part of the total Shapiro delay. As such it is 
entirely due to the mass of the companion star; in fact it 
is easy to verify that the term containing the mass of the 
pulsar emitting e.m. signals gives a contribution which 
remains constant during the orbital motion. Finally, the 
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third contribution is due to the gravito-magnetic field of 
the companion star. To lowest order, we obtain 



Ato ~ r\y + n sin i sin ( 



At, 



At M ~ 2A/ 2 In 



2J2 sin z cos t 



2r b 



n ■ r 



1 — n • f 



• 2 • • 2 

sin z sin 1 



(21) 



(22) 



(23) 



The contributions l|21l) and (|22|) to the time delay are 
in agreement with the standard results (see )> while 
equation H23|) gives the gravito-magnetic corrections. 

From the time delay ijL'211 we may extract the contribu- 
tion which varies during the orbital motion; if we use the 
equation of the orbit Q it can be written in the form 



At 



M 



-2Mo In 



1 — sini sin(u; + tp) 
1 + e cos if 

If we introduce the PK Shapiro parameters 

n = m 2) 

S = sinz, 



(24) 



(25) 
(26) 



the mass (or, gravito-electric) contribution to the time 
delay becomes 



At 



M 



-2ft In 



1 — S sin(o; + tp) 



1 + e cos ip 



(27) 



On the other hand, the gravito-magnetic contribution 
(|23[l changes continuously during the orbital motion, so 
that, similarly, we may write 



At* 



Ji sin z 
a(l-e 2 ) 



(cos(w + <p)) (1 + ecostp) 
1 — sinz sin(cj + (p) 



(28) 



If we introduce the PK parameter <S, and define a new 
PK parameter 



J = h 



(29) 



the gravito-magnetic contribution to the Shapiro time 
delay can be written in the form 



At*, = - 



JS 



a(l-e 2 ) 



(cos(w + >p)) (1 + ecosv?) 
1 — <Ssin(c<j + ip) 



(30) 



We notice that the new PK parameter J coincides with 
the intrinsic angular momentum of the source of the 
gravito-magnetic field. In particular, if we knew the rota- 
tion frequency of the source of the gravito-magnetic field 
(which is possible, if the latter is a visible pulsar, for in- 
stance), the new PK parameter could give information 
on its moment of inertia. 



The gravito-electric and gravito-magnetic contribu- 
tions to the time delay l|27ll and l|3U|) can be written in 
the form 



At 



At* 



A M F M (tp), 
AjFj(<p), 



(31) 
(32) 



where we have introduced the following constant " ampli- 
tudes" 



Aj 



-211, 



JS 



a(l-e 2 )' 



and the varying "phase" terms 



F M {<p) 



= In 



1 — 5 sin(o; + ip) 



1 + e cos ip 
(cos(w + </?)) (1 + e cos <p) 
1 — 5 sin(aj + tp) 



(33) 
(34) 



(35) 
(36) 



In particular we see that when S = sin?' = 1, Fj(tp) tends 
to diverge as far as u> + (p — > tt/2 (i.e. close to the con- 
junction position, when the impact parameter goes to 
zero). Fm too diverges in the same configuration and 
we see that the divergences have different "strength", 
the former being an inverse power, the latter logarith- 
mic. The reason for that difference is easily referable to 
the gravito-magnetic potential affecting Fj, which has a 
dipolar structure, and to the monopolar gravito-electric 
one affecting Fm- 

Even though the divergences have no physical mean- 
ing because the actual compact objects have finite di- 
mensions and the e.m. signals cannot pass through the 
center of the companion star, we see that the gravito- 
magnetic contribution is bigger for those systems that 
are seen nearly edge-on from the Earth, which are the 
ideal candidates for revealing the gravito-magnetic ef- 
fects. This is the case, for instance, of the binary sys- 
tem PSR J0737-3039, where, however, unfortunately the 
presence of a large magnetosheath zone makes the effec- 
tive impact parameter much bigger than the actual linear 
dimension of a neutron star [fj . That being said, in the 
following section we give numerical estimates for the con- 
stant amplitudes Am and Aj for the known binary pulsar 
systems. 



IV. NUMERICAL ESTIMATES 

The PK parameters related to the Shapiro time de- 
lay 1Z, S have been successfully measured with great ac- 
curacy in some binary pulsar systems, such as PSR B 
1913+16 (see ^2>EU)> an d the recently discovered sys- 
tem PSR J0737-3039 (see 0,0,(1). Indeed, the anal- 
ysis of these systems is very accurate, because of their 
favorable geometrical properties and, through the mea- 
surements of several PK parameters, they provided very 
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accurate tests of GR as confronted to alternative theories 
of gravity. In of , we studied the gravito-magnetic cor- 
rections to pulsar timing in a simplified situation, taking 
into account circular orbits only. Here, since we have 
generalized those results to arbitrary elliptic orbits, we 
may apply the formalism developed so far to all the bi- 
nary pulsar systems known up to this moment, in order 
to estimate the magnitude of the gravito-magnetic cor- 
rections to the time delay. In Table [I] Am and Aj are 
evaluated for the binary systems known until todavp9|. 
A few comments on how the table has been obtained: for 
those systems where the available data are not complete, 
we have chosen for the missing data the most favorable 
value (see the caption of the table). Furthermore, we 
have estimated the intrinsic angular momentum of the 
sources of the gravitational field supposing that the pro- 
genitor star was only a little bigger than the Sun, and 
that most of the angular momentum was preserved dur- 
ing the collapse, so that J 2 ~ J Q . 

From Table [I] it is clear that the gravito-magnetic con- 
tribution is much smaller than the mass contribution, as 
expected. However it is possible, at least in principle, 
to distinguish the former from the latter, on the basis of 
their different dependance from the geometric parame- 
ters. In fact, from l|2*l" |l .[|2*? )) and (j20J it is clear that the 
geometric and the gravito-electric contribution are sym- 
metric with respect to the conjunction and opposition 
points, while the gravito-magnetic contribution is anti- 
symmetric. As we pointed out in the previous paper |6j, 
if it were possible to identify conjunction and opposition 
points in the sequence of arriving pulses, this fact could 
be exploited for extracting the gravito-magnetic effect. 

Nowadays, the uncertainties in pulsar data timing are 
of the order of 10~ 6 ~ 10~ 7 s (see for instance |5|)- How- 
ever, as we pointed out above, gravito-magnetic effects 
can become larger if the geometry of the system is favor- 
able: in particular, when sini = 1, Fj tends to diverge 
close to the conjunction. On the other hand, we can 
give an estimate of the value of the geometrical parame- 
ters needed to make Aj of the order of magnitude of the 
present day uncertainties. So, if we assume J 2 — J©, in 
order to have Aj > 10 -7 s, we must have 



a (l_e 2 ) < a(1 , £2) <5x 10 4 m 



sin i 



(37) 



Hence, small orbits with great eccentricity allow, at least 
in principle, the measurement of the gravito-magnetic 
effects. It might be useful, also, to estimate the rate 
of decay of the orbit because of the emission of gravi- 
tational waves. If we assume, for the sake of simplicity, 
e = 0,Mi_= M 2 = 1.44M© and that a fulfills (E7J, we 
get (see [l^l) 



64 Mi Mi (Mi +M 2 ) 



< 2 x 10 7 s (38) 



Consequently, we may argue that these effects may be- 
come larger in the final phase of the evolution of the bi- 
nary systems, i.e. during their coalescence, even though, 
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TABLE I: Evaluation of the gravito-electric and gravito- 
magnetic contributions to the time delay. For the systems 
PSR J1756-2251RSR J1829+2456, PSR J1518+4904, PSR 
J1811-1736, PSR B2127+11C, since the present data do not 
provide the inclination of the orbit, we chose the most favor- 
able value for the calculations of Aj, i.e. sini = 1. Similarly, 
for the calculations of Am, we chose the best estimate for the 
mass of the companion star for the systems PSR J1829+2456, 
PSR J1518+4904, PSR J1811-1736, since the available data 
do not constrain it completely. 



in that phase, the weak field approach that we used in 
this paper would probably be rather poordemanding for 
different analysis techniques. 

The smallness of the gravito-magnetic delay would 
also require long data taking times so posing the prob- 
lem of the stability of the pulsar frequency. However a 
peculiarity which is not blurred by any drift or noise 
is the physical antisymmetry of the gravito-magnetic 
effect, which should emerge in the long period over all 
other phenomena. 



V. CONCLUSIONS 

In this paper, we have studied the effects of the gravi- 
tational interaction on the time delay of electromagnetic 
signals coming from a binary system composed of a radio- 
pulsar and another compact object. In particular, we 
have focused our attention on the gravito-magnetic cor- 
rections to the time delay due to the gravitational field 
of the binary system (Shapiro time delay). 

In doing so, we have generalized the results obtained 
in a previous paper, where we considered a simplified 
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situation, taking into account circular orbits only. Here 
arbitrary elliptic orbits are allowed. Furthermore, by fol- 
lowing a standard approach, we have expressed the time 
delay and its gravito-magnetic component in terms of Ke- 
plerian and post-Keplerian parameters. 

In particular, a new post-Keplerian parameter has 
been introduced, which coincides with the intrinsic an- 
gular momentum of the source of the gravitational field, 
and could give, in some cases, information on its moment 
of inertia. 

We have given numerical estimates of the amount of 
the gravito-magnetic corrections for all binary pulsar sys- 
tems known until today, and we have seen that, even 
though they are usually much smaller than the corre- 
sponding gravito-electric ones, they can become larger 
for those systems that are seen nearly edge-on from the 



Earth, which are the ideal candidates for revealing the 
gravito-magnetic effects in this context. 

Among the known binary pulsar systems, the one hav- 
ing the most favorable geometrical properties for the de- 
tection of the gravito-magnetic effect is PSR J0737-3039, 
which, unfortunately, has a large magnetosheath that 
keeps the magnitude of the gravito-magnetic correction 
below the detectability threshold. However, we cannot 
exclude, at the present discovery rate of new binary pul- 
sars, that other binary systems with favorable geometri- 
cal configurations can be found. This fact, together with 
the expected improvement of the sensitivity and precision 
of the timing of pulses, makes us cherish the hope that, 
in the future, it will be possible to measure the gravito- 
magnetic corrections to the time delay, and in particular 
the newly introduced post-Keplerian parameter. 
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